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ABSTRACT 

In this paper, considered a parabolic system of two quasi linear reaction-diffusion equations with the source 
and absorption task and the properties of self-similar sohitions of a system of quasi linear reaction-dijjiision equations 
for the source and absorption task. Self-similar system of equations is constructed by the method of nonlinear splitting. 
Estimates of the solutions and the free boundary that arises in this case are jound, which makes it possible to choose 
suitable initial approximations for each value of the numerical parameters. 
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INTRODUCTION 

Study of nonlinear mathematical models of various physical, biological, chemical and other phenomena 
and processes is one of the important directions of mathematical modeling. As examples, we note such physical 
theories as non-linear quantum mechanics, nonlinear electrodynamics and optics, non-linear theory of plasmas, 
nonlinear acoustics, nonlinear conduction, nonlinear diffusion, and other theories based on mathematical models 
which are nonlinear differential equations in partial derivatives. The study of the linear mathematical models of 
physical processes are easy to study, since the underlying linear differential equations developed general methods 
for their solution [1,2]. In applied tasks, the actual physical processes are nonlinear and for their adequate 
description, one should use nonlinear mathematical models. 

Nonlinear models of mathematical physics, describing phenomena and processes in a wider range of 
physical parameter changes and have a better capacity of information about these phenomena and processes. 
Linear models are typically special cases of nonlinear models. They can give only an approximate picture of the 
phenomenon under study without identifying the observed effects. Studies show that the nonlinearities change not 
only the quantitative characteristics of the processes, but the qualitative picture of their behavior. Interestingly, 
from the point of view of applications to study the following classes of nonlinear differential equations in which 
the unknown function and the derivative of this function consists of exponential way. Then, with the comparison 
theorems of solutions of this class can be extended. These types of nonlinearities are often encountered in 
problems of the theory of filtration, diffusion, thermal conductivity, magnetohydrodynamics, biological 
populations, the oil industry, etc. [3,4,5,6]. 
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These models more accurately describe the physics of the process and, therefore, their research shows that there 
are new effects related to the nonlinearity of the studied process. So was found the effects of the finite speed of propagation 
of perturbations [1], localization of solutions and different modes of processes. The first effect of finite speed of 
propagation of pcrturhations was obtained and applied to the problem of nonlinear thermal conductivity, in the work of 
Zeldovich and A. S. Kompaneetz[l], for the problem of nonlinear filtering in the work of G. I. Barenblatt[2,3], which are 
independently from each other, and got this effect. Further mathematical questions were studied in works of O. A. Oleinik, 
A. S. Kalashnikov [5,6,7], A. A. Samarskiy and his students [2,3-4], B. Knerr [3]. 

Statement of the Problem 


In the field Q=[(t,x): 0<t<°°, xeR} considers a parabolic system of two quasilinear equations of source and 
absorption 


3mj d | 
dt dx 
du 2 3 
dt dx 


D x u°' 

D 2 u 2 2 + k,(t)u, • (l- uf -), 

dx J 


( 1 ) 


«i| f =o-«io(*). U 2 \ t=0 , (2) 

Which describes the process of heat distribution in a nonlinear two-component environment, the diffusion 
coefficient of which is equal to D X U°' and D 2 U 2 2 , &\,<J 2 ,f}\,f} 2 - positive real numbers, U x = U t (t,x) > 0, 
U 2 = U 2 (t,X ) > 0 - desired solution. 


Estimates of solutions are found, which makes it possible to choose suitable initial approximations for each value 
of numerical parameters. 

Construction of a Self-Similar System of Equations 

Self-similar system of equations is constructed by the method of nonlinear splitting[4]. 


— (f°' + u ( f f f&) = 0 

—( fT- %++ Mh ~ yi /i A ) = o. 

d£ 2 d£ 2 d% Jl yiJ ' 


where 


Mi = 


1 

(i -m)’ 


To construct the upper solution of the system (1) - (2), we introduce functions 

/(#) = A(a - # 2 ) + ">, f 2 (£) = B(a - ? ) + " 2 , 


where n, =-, n 2 =- , (b) + = max(0, b ), 

<j\ <J 2 
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~^k i (q)dq _ 

U\ (t, x) < u 1+ {t, x) = e 0 / (#), 

-jk 2 (q)dq _ 

u 2 {t,x)< u 2+ (t,x) = e 0 f 2 (£), 


/,(#). /(#)« 


*■(<) 


■the functions defined above. 


Cross-Diffusion System of Two Quasilinear Reaction-Diffusion Equations of the Source and Absorption 

In the domain Q={(t,x): 0<t<°o, xeR} considered a parabolic system of two quasilinear reaction-diffusion 
equations for the source and absorption 


du { 

~dt~ 

du 2 

lk 


d_ 

dx 

d_ 

dx 


D u a ' 

U \ U 2 

OX y 

D 2 M I 

ox 


+ 


+ 


k { (t)u { [\-u 2 ), 
k 2 (t)u 2 (l-uf 2 ), 


(3) 


M i| f =o- M io(*)> ^2(1=0 u 2 y(x ), 

which describes the process of process of absorption in a nonlinear two-component environment, diffusion 
coefficient whose are D j uX' and D 2 U° 2 , <7 { ,(7 2 ,/3 { , j3 2 ~ positive real numbers, W, = U { (t,x) > 0 , 
U 2 = U 2 (t,x) > 0 - sought solutions. 


Let 7 2 £7, > 1, 7,£7 2 = f,£7,, £7 2 (b 2 +1) + /3 2 (b { + 1) = (7, +1) + /?, (fc 2 + 1) , C,- > 0 . In this case, 

assuming in (3) 

W, (T(t),x) = f .,(#), ^ =1 JCI /(r/ 72 , «=7,2, 


and taking into account that the equation for W- (X, +) without lower terms always have a self-similar solution in 
the case 1 — /j_, £7 + 0, obtain the system 




(4) 


where 


A = 


1 

1 - ' 


To construct the upper solution of the system (4), we introduce functions 
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f(f = A(a- f )/, /, (f = B(a - f ! )f , 


where 


n, 


(<7, +2) 


l-(cr,+l)((J 2 +l)' ' ?2 1-(<7 j +l)(cr 2 + 1) 


( ^ + 2) _, (b) + = max(0,Z?). 


and 


f 2 ^ = -2B^' r2 tfeC(0~), 
dg 

ff 2 & = -2A^ 2+l r ^f 2 GC(0, oo) 




d_\ 7+ dfi 

d %l 2 


d 


Ja 2 df 

d^{ 1 


= -2^ 


= -2ftA° 


C ,+1 


/i + £ 




d% 


A and B are chosen from the system of nonlinear algebraic equations 


r2 B a ' +l = 1 / 2 , 


ri A^ 2+I =1/2, 


and 


/?, =1 / n 2 , (3 2 = 1 / n,. 

Then in the domain Q according to the principle of comparing solutions we have 

t 

-^k^drj _ 

u x (t,x)<u l+ (t,x) = e 0 f x (£), 

t 

-^k,(rj)dri _ 

u 2 (t,x)<u 2+ (t,x) = e 0 / 2 (f), 

where /,(#>, /+a w m -previously defined function. 

Slow diffusion.Case /7, > 0, U 2 > 0, n > 0 (slow diffusion). Using the method of nonlinear splitting for 
the solution of the equation (4) obtained the following functions 

e^) = (a-f-):\ 0 2 ({) = (a-f-):\ 

where 
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(<T, + 2) (<T 7 + 2) 

=— 1 --■ =— --, n = l-^Cj + l)(cr 2 + l), 


n 


n 


a > 0, ( v) + = max ( y, 0) , ^ < d . It is known that for the global existence of a solution of the task (4) of the 
functions f\ (<f ) must satisfy the following inequality: 




and 




A = 1 / n 2 , A = 1 /«i. 

It is shown that the functions ^ (^), (^) will be asymptotic finite solutions. 

TeopeMa 1. Finite solution of the problem (4) at ^ —> Cl has asymptotics ( f ) , i — 1,2 . 

Fast diffusion.Case Wj > 0, n 2 > 0, n < 0 (fast diffusion). For (4) there is 

Z(f) = (a + f) \ Z 2 (£) = ( a + f~) ,l \ 


where 


-ta+2) 


n.=^^.n,= ~ (<T - + 2) . k = (c r,+l)(<r 2 +l)-l. 


a >0. 

Theorem 2. At C —> +°° the vanishing at infinity solution of the problem (4) has asymptotic f\(Z)~Z, ( < f) ■ 

And also in the domain Q={(t,x); 0<t<°°, xeR} considered parabolic system of two quasilinear equations of the 
source and absorption task 


du x 

_ 3 " 

d t 

dx 

du 2 

_ d " 

dt 

dx 


(a u u'f +a n u’f)^ + (b u uf +VD^ 
ax ~ ax 


+ k x (t)u x (\-U2 ) 


(a 21 u x + a 22 u'f)-\- + (b n u x +b 22 u 2 + k 2 (t)u 2 (l-uf 2 ) 
ax ~ ~~ ox 


l -, b X j - positive real numbers, p x ,p 2 > 0, U x = U x (t , x) > 0, U 7 = U-,(t, x) > 0 - desired solution. At 


(5) 


a u* 0’ =0 or a X j — 0 , bjj + 0 mathematical model (5) is a system of the type reaction-diffusion with diffusion 
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coefficients a U m > 0, b U m > 0 . In the case where at least one of the coefficients a= + 0 and b sl + 0 (the sign can 

y i y 1 y y 

be arbitrary), system (5) is cross-diffusion (mutual-diffusion for i,j=l,2). 

Qualitative properties of the considered problem are studied by constructing self-similar system of equations for 

(4). 


Then, the obtained self-similar 

d 


system 


d 


df\ 

i tU f m i u f m \ df 2 




(«n/r + a l 2 f 2 m )-^ + (b l J 1 m + h nf?)j^ 

(a,J~ +a j-) l f+(bj: +bJ-)% 
dg dg 


+~+w, f ,a- f f)=o. 

+~+v,m-ff'-)=o. 


(6) 


and r, = r, (t) is chosen as 


(T + T)- m+1 

—y ] m +1 


if -y x m + 1^0, 




ln (T+r), 


(T + t), 


if - y x m + 1 = 0, 
if m = 0, 


if y 2 m x = y,m 2 . 

System (6) has an approximate solution of the form 

/i = A(a-bf~) + m , f 2 =B(a-bf-)f 2 (y) + = max(0, y). 


Case I 


y, =y 2 > V) 


m 


and the coefficients A and B are determined from the solution of the systems of nonlinear algebraic equations 

j[a n A" + a n B m \l),A I [/.,,A" +b u B m \t] 1 B■ (-2b) = A. 
lk,A" + a n B"\ T),A + [ b n A " + b 22 B m \f) 2 B ■ (-2 b) = B. 


Case II 


a n = 0; a l2 + 0; b u = 0; b l2 = 0; 
a 21 = 0; a 22 = 0; b 2[ + 0; b 22 = 0; 

1 1 

>72=-> %=—• 
m m 

And the coefficients A and B are determined from the solution of the system of nonlinear algebraic equations 
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a u B m ri x (-2b) = \. 


b 2l A n -ri 2 {-2b)=\. 


Case III 


a n =0; a l2 =0; 
a 2l ± 0; a 22 = 0; 


b u = 0; b n ^O; 

^21 — ^22 — 


m + 2 _ m + 2 

( m + 1) 2 -1 _ ( m + 1) 2 -1 


And the coefficients A and B are determined from the solution of the system of nonlinear algebraic equations 

\b u B m+l -7l 2 {-2b) = A, 

\ a 2l A m+1 ■ rj x (~2b) = B. 

Types of Self-Similar Regimes with Peaking 

There can be three types of self-similar modes with peaking: NS, S, and LS. At 0 < /? < , i = 1,2 HS mode 

is realized. Studies have shown that a self-similar problem in this case has a unique eigenfunction that decreases 
monotonically on a segment with a maximum at the center of symmetry. A self-similar solution is a wave whose amplitude 
and front increase in the blast mode. 


At /? ; = ITlj, i = l,2takes place S-mode. Self-similar solution is a non-stationary dissipative localized structure. 

Inside the area of localization, the number of individuals increases in regime with aggravation, and beyond remains equal 
to zero. 

Self-similar solution in LS-mode, is a non-stationary dissipative structure; all points of which move to the center 
of symmetry, the solution at T = — T turns into infinity in a single point - the center of symmetry. Self-similar solutions 
can exist when /3j > in t , i = 1,2. 


The following are the results of numerical experiments for different values of the parameters. 

To determine the number of eigenfunctions (EF) self-similar problem in LS - mode, and studies of the nature 
of their dependence on the parameter was conducted bifurcation analysis of the solutions. The first ten eigenfunctions 

continued on the parameter /? ( for a fixed parameter m t ( m { =2). 


www.tjprc.ors 


editor@tjprc. org 





562 


D. K. Muhamediyeva 



Figure 1: Results of Computational Experiment 


CONCLUSIONS 

The analysis showed that the eigenfunction with the number j — 2,3,4,... exists in the interval 
/7t ( . < < /3-, where 
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Pij = + m t /(j - 1), i = 1,2, j = 2,3,4,... (7) 

Values /J = m t h A = Pu are the bifurcation points in which the EF ceases to exist. The first EF exists for 

any value j3 t > m { , i — 1,2 . It follows from (7) that when > /3 i2 = 2m t , i — 1,2 self-similar task in LS mode can 

have only one eigen function. 

The greater the number of SF, the narrower the interval of the parameter /3 t in which it exists. 

The number of eigenfunctions N that has a self-similar problem for the given f3 t and iYl t is determined by the 
formula, where Cl = (3 ; K(3 t — Hi t ) , i = 1,2 : 


J[a i ],(/’ a is not an integer, 

[a,. - 1, if a is integer. 

The sequence of bifurcation points that define the right boundary of the domain of existence of eigenfunctions 
is an infinite sequence that converges to a point 1 + = oo op, which is the common left boundary of the intervals of 
existence of all eigenfunctions in the LS-mode. 
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